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ABSTRACT 

Recent theoretical progress in the evaluation of Higgs boson production at the 
■ LHC within the Standard Model and its minimal supersymmetric extension is 

Q\ \ reviewed. In particular the two-loop QCD corrections to the squark loop con- 

tributions to scalar Higgs production in the MSSM and soft gluon resummation 
^ | effects in Standard and SUSY Higgs production via the gluon fusion mechanism 

are discussed. 



qq ■ 1. Introduction 

The search for Higgs particles lil is one of the most important endeavours for 
; future high energy e + e~ and hadron collider experiments. The Higgs boson is the 

only particle of the Standard Model (SM) which has not been discovered so far. The 
£T) ■ direct search at the LEP1 experiments via the process e + e~ — ► Z*H yields a lower 

bound on the Higgs mass of ~ 65.2 GeV i. Theoretical consistency restricts the Higgs 
mass to be less than ~ 700 GeV i. The dominant Higgs production mechanism at 
the LHC, a pp collider with a cm. energy of 14 TeV, is the gluon fusion process 
gg — > H, which is mediated by a heavy quark triangle loop at lowest order @. As an 
important step to increase the theoretical precision, the two-loop QCD corrections 
have been calculated, resulting in a significant increase of the predicted total cross 
section by about 50-100% H. The dependence on the unphysical renormalization and 
factorization scales decreased considerably by including these next-to-leading-order 
^ . (NLO) corrections, resulting in an estimate of about 15% for the remaining scale 

sensitivity i. It is important to note that the NLO corrections are dominated by soft 
gluon radiation effects. 

The minimal supersymmetric extension (MSSM) is among the most attractive 
extensions of the SM. It requires the introduction of two Higgs doublets leading to 
the existence of five scalar Higgs particles: two scalar CP-even h, H, one pseudoscalar 
CP-odd A and two charged bosons H ± . This Higgs sector can be described by 
two parameters, which are usually chosen to be tg/3, the ratio of the two vacuum 
expectation values, and the pseudoscalar Higgs mass Ma- Including higher-order 
corrections to the Higgs masses and couplings up to the two-loop level, the mass of 
the lightest scalar Higgs particle h is restricted to be smaller than ~ 130 GeV @. 
The direct search at LEP1 sets lower bounds of about 45 GeV on the masses of the 
MSSM Higgs bosons i. The dominant neutral Higgs production mechanisms at the 
LHC are the gluon fusion processes gg — > h,H,A and Higgs-strahlung off bottom 
quarks, gg, qq — > bbh, bbH, bbA, which become important only for large tg/3 i. The 



1 



coupling of the neutral Higgs particles to gluons is again mediated by top and bottom 
loops, with the latter providing the dominant contribution for large tg/3, and squark 
loops, if their masses are smaller than about 400 GeV i. The two-loop (NLO) QCD 
corrections to the quark loop contributions to the gluon fusion mechanism have also 
been calculated @i, and conclusions completely analogous to the SM case emerge. Soft 
gluon radiation effects again provide the dominant contribution to these corrections, 
for small tg/3. 



2. Gluon Fusion: Squark Loops 



Recently the QCD corrections to the squark loop contributions to the cross sec- 
tions a(pp — > Ti. + X) of the fusion processes for the neutral CP-even Higgs particles 
H = h,H 

99 -> U{g) and gq -> Hq, qq -> Kg (1) 

have been calculated. Because of CP invariance, squark loops do not contribute to the 
production of the CP-odd Higgs boson in lowest order. The QCD corrections from 
squark loops were evaluated in the heavy squark limit, where the calculation can 
be simplified by extending the lowest-order low-energy theorems to two loops. 
This limit should be a good approximation i for the production of Higgs particles 
with masses smaller than twice the squark masses. For simplicity, we will restrict 
ourselves to the case of degenerate squarks where mixing effects are absent. Also in 
this case, scalar squarks do not contribute to the production of the CP-odd Higgs 
boson A at NLO. 

The LO cross sections for CP-even Higgs production at the LHC are given by 

dC" 



CT L0 (PP — n + X) 



dsn 



(2) 



H 



with dC" jdr-n the gluon luminosity at r-n = M^/s and s the total cm. energy. The 
parton cross sections are built up from heavy quark and squark amplitudes, 
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The normalized scalar quark and squark couplings to the CP-even Higgs bosons, gV; 
can be found in Ref.H. The sums run over t, b quarks and the left- and right-handed 
squarks Ql, Qr, which in the absence of mixing are identical to the mass eigenstates. 
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The QCD corrections to the gluon fusion process, eq. ([]]), consist of virtual two- 
loop corrections and one-loop real corrections due to gluon radiation, as well as con- 
tributions from quark-gluon initial states and quark-antiquark annihilation. The 
renormalization program has been carried out in the MS scheme for the strong cou- 
pling constant and the parton densities, while the quark and squark masses are defined 
at the poles of their respective propagators. The result for the cross sections can be 
cast into the form 



H 



1 + C H (t Q ,Tq)^ 



dC" 

TH^ + ^ + AaK + Aa*. (5) 



The coefficient denotes the virtual two-loop corrections regularized by the infrared 
singularities of the real gluon emission. The terms Aajj = g, q) denote the finite 
parts of the real corrections due to gluon radiation and the gq and initial states. 
The expressions for the t, b quark contribution can be found in Refs.Era. 

The calculation of the QCD corrections has been performed by extending the 
low-energy theorems to scalar squarks at the two-loop level. For a light CP- 
even Higgs boson, these theorems relate the matrix elements of the quark and squark 
contributions to the Higgs-gluon vertex to the gluon two-point function. In the 
following we consider only the pure gluon exchange contributions, which are expected 
to be the dominant ones; for heavy enough gluinos, the two-loop corrections due to 
gluino exchange should be small, since they are suppressed by inverse powers of the 
gluino mass. The final result for the squark contributions to the 7i coupling to gluons 
can be expressed in terms of the effective Lagrangian [v = (v^Gp) -1 / 2 ] 



96 PoM/a. H ^ u <*s H L 25 a 



1+ 

O 7T 



(6) 



where we have used the anomalous squark mass dimension 0, 7g = 4ct s / (37r), and the 
squark contribution to the QCD (3 function 0, /3g(a s ) = a 2 s /{\2n) [1 + lla s /(27r)]. 
Starting from the Lagrangian of eq. @, the effective QCD corrections due to real 
gluon emission and the gq/qq initial states have to be added. These corrections are 
identical to the corresponding corrections to quark loops @ in the heavy quark limit. 

The QCD-corrected squark loop amplitudes have to be added coherently to the 
corrected t, b loop amplitudes, whose full mass dependence is known. To obtain a 
more reliable prediction for the total cross sections, the resulting amplitudes for the 
squark contributions have been normalized to the lowest-order amplitude in the limit 
of large squark masses. These ratios are then multiplied by the lowest-order amplitude 
including the full squark mass dependence. The heavy squark limit is then expected 
to be a very good approximation for Higgs masses below the QQ* threshold, as in 
the corresponding case of top quark contributions 1. The final results for the partome 
cross sections defined in eq. (W) can be found in I. 
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Fig. 1. Ratio of the QCD- corrected cross section to the lowest order result for o~(pp —fTL + X) 
for tg(3 — 1.5 and 30. Solid lines include t,b and squark contributions (with wlq — 200 GeV); 
dashed lines include only the t,b contributions. We have taken m,b ^,5 GeV, mt — 176 GeV and 
a s (Mz) = 0.118. We use next-to-leading order GRV parton densities^ and take the renormalization 
scale [i and the factorization scale M equal to M h / H . 



In Fig. 1, we present the i^-factors for the QCD corrections to the production 
of the CP-even MSSM Higgs bosons as functions of the Ti masses for the LHC at a 
cm. energy \fs = 14 TeV with (solid lines) and without (dashed lines) the squark 
contributions. The ^-factors are defined as the ratios of the QCD-corrected and 
lowest-order cross sections, using next-to-leading order a s and parton densities in 
both terms. A common value itlq = 200 GeV has been used for the left- and right- 
handed squark masses. This value is identified with the SUSY scale of the MSSM 
couplings and Higgs masses. The QCD corrections to the squark loops are large, 
approximately of the same size as the QCD corrections to the quark loops. They 
enhance the cross sections by a factor between 1.6 and 2.8. However, if the lowest- 
order cross sections are convoluted with lowest-order a s and parton densities, the 
.fT-factors are reduced to a level between 1 and 2. It can be inferred from Fig. 1 that 
the inclusion of squark loops in the production of both CP-even Higgs particles h and 
H does not substantially modify the ii'-factors compared to the case where squark 
loops are absent. Thus, to a good approximation, the effect of the squark loops in the 
gluon fusion mechanism is quantitatively determined by the lowest-order cross section 
(including squark loop contributions), multiplied by the known i^-factors when only 
the t, b quark contributions i are taken into account. 

3. Soft Gluon Resummation 

As a next step a resummation of soft gluon radiation effects in Higgs boson pro- 
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duction via gluon fusion has been performed. If the i^-factor, evaluated in the heavy 
top quark limit, is multiplied by the full massive lowest-order cross section, the result 
approximates the SM gluon-fusion cross section at NLO within 10% and the MSSM 
cross sections within 25% for tg/3 < 5 0. In the heavy quark limit the partonic cross 
sections factorize as 

with the coefficient ctq defined in eq. 



and 



0n 



12872 



7T 



Y,9QT Q f(r Q ) 



where gq{4> = h, H, A) denote the modified top Yukawa couplings normalized to the 
SM coupling, which are given in i. In the following we will neglect the gq, qq initial 
states, which contribute up to ~ 10% at NLO, and the squark loops, by choosing 
the common squark mass as itlq = 1 TeV. The factor in eq. fl7|) stems from the 
effective coupling of the Higgs boson to gluons in the heavy top quark limit, which 
can be obtained by means of low-energy theorems They are given" by 
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11 a s (m 2 t ) 3830 - 201 n/ 
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TV 



+ 0{al) 



(9) 
(10) 



where a s is the strong coupling constant in the MS scheme including n/ = 5 flavours in 
the evolution, i.e. the top quark is decoupled. The scale of the higher-order corrections 
to this effective coupling has to be identified with the top quark mass m t 

We will now construct a resummed expression for the correction factors p^ by 
means of the methods described in Ref. 0. Near the elastic edge of phase space the 
Higgs cross section in the infinite mass limit may be factorized into hard, soft and 
jet functions, in complete analogy with the Drell-Yan cross section. Following the 
arguments of Ref. Ill this leads to the Sudakov evolution equation: 




M l f 2 
v ' p 2 



The solution to the moments^ of eq. (0) is given by 



' M 2 
N } -jt } a s (p 2 ) 



exp 



(12) 



a The factors kh,h include the top quark contribution at vanishing momentum transfet- which differs 
from the top quark contribution to the MS (3 function by a finite amount at 0(a^) EH. 
& We define the Mellin-moments in the usual way as f(N) = L dzz N ~ 1 f(z). 
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which may be expanded perturbatively. In this way the one-loop coefficients of 
the evolution kernel can be straightforwardly determined from the explicit NLO 
calculation 0. In the following we will define three approximations: scheme a includes 
only the leading logiV contributions, scheme (3 all terms of O(N ), while scheme 7 
contains terms of O (log N/N) in addition, which arise from collinear gluon-gluon 
splitting and are thus universal. [However, they are not covered by the present level 
of factorization theorems 0, which leads to eq. (|TT|).] 

In order to obtain a finite expression for the correction factors we now have to 
renormalize the strong coupling constant and perform mass factorization of the gluon 
distributions. Both objects have been defined in the MS scheme. The final results 
for the moments of the renormalized correction factors can be found in @. They have 
been expanded perturbatively up to NLO and NNLO, thus yielding a quantitative 
estimate of the unknown NNLO corrections to the gluon-fusion mechanism 0. 
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Fig. 2. a) Exact and approximate two- and three-loop correction factors, convoluted with the NLO 
gluon-gluon luminosity dC a ^ LO / dr , in the heavy top-mass limit. The results for the three different 
schemes are presented as a function of the scalar Higgs mass Mh, using NLO CTEQ4M parton den- 
sities llil and a s [A&L = 202 MeV]. b) Hadronic NLO K- factor using LO CTEQ4-L parton densities 

and a s [A ( H = 181 MeV] for the LO cross section and including the NLO contributions from Kh- 

Fig. [|a shows the NLO and NNLO expansions of the scalar correction factor, con- 
voluted with NLO gluon densities, as a function of the Higgs mass. A similar picture 
emerges for the pseudoscalar Higgs boson. The NLO correction factor in scheme 7 
coincides with the exact NLO result within less than 5%, while schemes a and (3 
fail to approximate the exact NLO calculation reliably. From an analogous analysis 
of the Drell-Yan process we get strong confidence that scheme 7 also approximates 
the NNLO expansion with reliable accuracy 0. Thus it might be expected that the 
curve labelled 72 yields a reasonable approximation of the unknown NNLO correc- 
tions to the gluon-fusion process. The NNLO correction factor amounts to 2.7-3.5 in 
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the relevant Higgs mass range, which is very large. However, the physical f^-factor 
requires the LO result to be convoluted with LO parton densities and strong coupling 
and the NNLO expansion with NNLO quantities. This leads to a strong reduction 
of the NLO i-T-factor in comparison to the naive correction factor, as can be inferred 
from Fig. |2|b. Due to the lack of NNLO parton densities no physical prediction of the 
Higgs production cross section at NNLO can be made. From the situation at NLO 
one might expect a significant decrease of the NNLO K-i actor compared to the result 
of Fig. ||a. 

An alternative way to reach a consistent prediction of the gluon-fusion cross sec- 
tion at NNLO may be provided by extracting the gluon luminosity from another 
production cross section at the LHC, e.g. gg — > ti, the resummed expression of which 
also has to be expanded up to NNLO. This result may then be inserted in the convo- 
lution integral for the Higgs production cross section. This method would thus allow 
us to obtain approximate NNLO parton densities for physical predictions of produc- 
tion cross sections at the LHC. In the same way we could also perform predictions of 
the full resummed results. This, however, requires a prescription, how to regularize 
the infrared renormalon singularity El, which appears in the integrals of the running 
strong coupling constant. 
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Fig. 3. Scale dependence of the Hiaas p roduction cross section for two values of the Higgs mass 
Mh. NLO CTEQ^M parton densities Ej and strong coupling /A^- = 202 MeV] have been used in 
all expressions, so that the NNLO results do not correspond to the physical NNLO cross sections. 

Fig. H presents the scalar Higgs production cross section at LO, NLO and NNLO 
as a function of the renormalization/factorization scale in units of the Higgs mass 
for Mh = 500 GeV. Again, all orders of the cross section have been evaluated with 
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NLO parton densities and strong coupling, and thus the NNLO curves do not cor- 
respond to the physical cross section at NNLO. The dashed lines present the NLO 
and NNLO results of scheme 7, where the NNLO contribution has been added to the 
exact NLO result. The solid NLO curve corresponds to the exact NLO cross section 
in the heavy top quark limit, while the solid NNLO line includes the exact scale de- 
pendence at NNLO, which has been obtained from the exact NLO result by means of 
renormalization group methods @. The full NNLO curve has been normalized to the 
72 curve at fi — M — Mh- Fig. ^| clearly indicates a strong stabilization of the scale 
dependence at NNLO, which develops a broad maximum around the natural scale 
fi = M ~ Mh in contrast with the NLO curve, which is a monotonic function of the 
scales. Thus with NNLO parton densities it would be possible to predict the Higgs 
production cross section at the LHC accurately with small theoretical uncertainties 
from the remaining scale dependence. However, the non-perturbative uncertainties 
related to the renormalon singularity of the full resummed result may be sizeable. 
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